Maximal air bubble entrainment at liquid drop impact 
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At impact of a liquid drop on a solid surface an air bubble can be entrapped. Here we show 
that two competing effects minimize the (relative) size of this entrained air bubble: For large drop 
impact velocity and large droplets the inertia of the liquid flattens the entrained bubble, whereas for 
small impact velocity and small droplets capillary forces minimize the entrained bubble. However, 
we demonstrate experimentally, theoretically, and numerically that in between there is an optimum, 
leading to maximal air bubble entrapment. Our results have a strong bearing on various applications 
in printing technology, microelectronics, immersion lithography, diagnostics, or agriculture. 
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The impact of liquid droplets on surfaces is om- 
nipresent in nature and technology, ranging from falling 
raindrops to applications in agriculture and inkjet print- 
ing. The crucial question often is: How well does the 
liquid wet a surface? The traditional view is that it is 
the surface tension which gives a quantitative answer. 
However, it has been shown recently that an air bubble 
can be entrapped under a liquid drop as it impacts on the 
surface Also Xu et al. 0, §] revealed the impor- 

tant role of the surrounding air on the impact dynamics, 
including a possible splash formation. The mechanism 
works as follows The air between the falling drop 

and the surface is strongly squeezed, leading to a pres- 
sure buildup in the air under the drop. The enhanced 
pressure results in a dimple formation in the droplet and 
eventually to the entrapment of an air bubble (figured^,). 
The very simple question we ask and answer in this pa- 
per is: For which impact velocity is the entrapped bubble 
maximal? 

Our experimental setup is shown in figure [Tp and is 
similar to that of ref. @ where it is described in detail. 
An ethanol drop impacts on a smooth glass surface after 
detaching from a needle, or for velocities smaller than 
0.32 m/s, after moving the needle downwards using a lin- 
ear translation stage. A high-speed side view recording is 
used to measure the drop diameter and velocity. A syn- 
chronized bottom view recording by a high-speed color 
camera is used to measure the deformed shape of the 
liquid drop. Colored interference patterns are created 
by high-intensity coaxial white light, which reflects from 
both the glass surface and the bottom of the droplet. 
Using a color-matching approach in combination with 
known reference surfaces, the complete air thickness pro- 
file can be extracted (shown in figure[TJ:). For experi- 
ments done at larger impact velocities (U > 0.76 m/s), 
we use a pulse of diffused laser light triggered by an op- 
tical switch. The thickness of the air film at the rim is 
assumed to be zero, and the complete air thickness pro- 
file can then be obtained from the monochromatic fringe 
pattern. From these measurements we can determine the 
dimple height, Hd, and the volume of the entrained bub- 
ble, Vb, at the very moment of impact. This moment 



is defined by the first wetting of the surface, i.e., the 
moment when the concentric symmetry of the interfer- 
ence rings is lost. To calculate the bubble volume Vb, 
we integrate the thickness profile of the air layer trapped 
beneath the drop. Alternatively, we can also measure 
the volume of the trapped bubble after impact when the 
liquid already wets the surface. Both measurements pro- 
vide the same results. In the present article, we use the 
first approach. 

The results are shown in figure [H Clearly, both dim- 
ple height at impact and the size of the entrained bub- 
ble have a pronounced maximum as function of the im- 
pact velocity U. The corresponding impact velocity for 
which the air entrainment is optimal is U — 0.25 m/s 
for an ethanol droplet of radius R = 0.9 mm (or the 
Stokes number St = 1 x 10~ 4 ). While length scales 
are given in multiples of the droplet radius R, follow- 
ing Brenner and coworkers 3, 6] we express the impact 
velocity U in terms of the (inverse) Stokes number St, 
defined with the dynamic air viscosity r\ g and the liq- 
uid density pi as St = r]g/(piRU) = p g jpiRe~ x , where 
Re — PgRU /rj g is the standard Reynolds number. A 
further relevant parameter of the system is the surface 
tension 7, which can be expressed in terms of the We- 
ber number We = piRU 2 /j or in terms of the capillary 
number Ca = r/gU/j = St ■ We. 

We compare and supplement our experimental findings 
on the dimple height at impact and the entrained bubble 
size to numerical results. The numerics consists of an 
axisymmetric boundary integral (BI) simulation for the 
liquid droplet (i.e., the droplet is assumed to obey poten- 
tial flow) , coupled to a lubrication approximation of the 
Stokes equation 
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that describes the viscous, incompressible gas flow under 
the droplet @,[lMl- Here, P g (r,t) is the gas pressure, 
while u r is the radially outward velocity in the gas par- 
allel to the surface (figure QJi). Note that the gas flow 
under the droplet is indeed viscous: An upper bound for 
the Reynolds number relevant for the lubrication flow 
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FIG. 1: Experimental characterization of air bubble entrapment, (a) Sketch of dimple formation (not drawn to scale) just prior 
to impact, (b) Schematic of the experimental setup used to study droplet impact on smooth surfaces. An ethanol droplet of 
typical radius R = 0.9 mm falls on a glass slide of average roughness 10 nm. The impact velocity is varied by varying the falling 
height of the droplet. For very small velocities below 0.31 m/s, the droplet is fixed at the tip of 0.4 mm-diameter capillary that 
is vertically translated downwards at a constant velocity. The bottom view is captured by a high-speed color camera (SA2, 
Photron Inc.). The camera is connected to a long working-distance microscope and a 5x objective to obtain a 2mm field of 
view, (c) An example of an interference pattern and the extracted air thickness profile. Note the difference in horizontal and 
vertical length scales. 



gives UHd/vg ~ 0.1 for the highest impact velocity, and 
is typically much smaller for most of our experiments. 

We now give more details on the numerical simula- 
tion: The velocity field inside the droplet is described 
with a scalar velocity potential <j>, obeying the Laplace 
equation V 2 </> = 0. The axisymmetric droplet contour is 
described using cylindrical coordinates r, z and is solved 
numerically by using the BI method; the simulations are 
based on the numerical code described by [T3 - [l6j . This 
BI simulation is an alternative way of solving the sys- 
tem of equations, compared to the method applied by 
10], in which case a Hilbert transform method was ap- 
plied. In contrast to [HI, we do not solve the complete 
Navier-Stokes equations, but do include dynamics of the 
air layer below the drop. The dynamic boundary con- 
dition on the droplet contours is given by the unsteady 
Bernoulli equation, 
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Here t is time, g the acceleration of gravity, z the abso- 
lute height, n(r,t) the interface curvature, and the 
far-field pressure. The key dynamical quantities in ([2]) 
are the gas pressure P g (r, t) and the interface curvature 
K,(r,t). The curvature is related to the dimple profile 
H(r,t) by the geometric relation 
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To close the problem, an additional equation is provided 
by the lubrication approximation for the viscous gas flow 
at the bottom of the droplet, 



dH(r,t) _ ld_ 
dt r dr 
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with boundary condition P g \ r= R = Poo] the gas pressure 
at the top of the droplet is set to atmospheric. Contrarily 
to [Io|, we do not incorporate effects of compressibility 
of the gas, since, following the analysis of |l2j, there is 
little influence of compressibility in the regime that is 
studied here. The initial conditions for the simulations 
consist of a spherical droplet with radius R with a down- 
ward velocity U . The initial height is taken sufficiently 
high for the viscous lubrication pressure to be still neg- 
ligible (~ 10 /mi). The number of nodes on the droplet 
surface for which the BI equations are solved is of order 
100, with node density increasing for r — > 0. The num- 
ber of nodes and the size of the time steps vary during 
the simulation, as a function of the local gap height and 
velocity of the droplet contour. The size of a time step 
is of order 10 ns. For any number of nodes, the coupling 
between gap height and pressure profile breaks down for 
some small value of H, since the pressure diverges at 
vanishing thickness of the air layer. Consistent with the 
experimental resolution we continue our simulations un- 
til the minimum gap thickness reaches 0.4 /tm, while en- 
suring that our algorithm remains accurate. This is the 
moment at which the values for Hd and V& are extracted. 

The results of the numerical calculations are shown 
in figure [21 together with the experimental data, show- 
ing very good agreement with our experimental results: 
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in particular, we observe the pronounced maxima in 
the dimple size and in the entrained bubble volume at 
the optimal Stokes number St a . In the numerically ob- 
tained bubble volume, we observe a jump exactly at the 
crossover regime. This jump originates from a change in 
the shape of the dimple. Figure [3] compares the experi- 
mental and numerical dimple profiles for an impact veloc- 
ity at the crossover regime (C/ =0.2 m/s) and an impact 
velocity in the inertial regime (U—0.7 m/s). While the 
profiles are in excellent agreement in the inertial regime 
(both volume and dimple height), the numerical profile 
develops a "double dimple" at the lower impact speed. 
This variation in shape results in the jump observed for 
the numerical bubble volumes in the crossover regime 
(see figure EJd). In all cases, however, the dimple height 
Hd is in quantitative agreement without any adjustable 
parameters. 

Numerical and experimental results together suggest 
scaling laws H d /R ~ St 2 ' 3 and V b /R 3 ~ St 4 ? 3 for smaller 
Stokes numbers, while H d /R ~ St- 1 / 2 and V b /R 3 ~ 
St- 1 for larger Stokes numbers. We will now theoret- 
ically derive these scaling laws. For small St we follow 
and extend ref. @, [ll|: The horizontal length scale L 
of the dimple extension (see figure [T^) follows from ge- 
ometrical arguments as I ~ ^/HdR, and u r from mass 
conservation as u r ~ UL/Hd- The Stokes equation ^ 
suggests P g ~ Lr/giir/Hj as estimate for the gas pressure 
below the falling drop at touch-down. The liquid pressure 
Pi can be estimated from the unsteady Bernoulli equa- 
tion: dimensional analysis for the deceleration timescale 
Hd/U and for the potential in the liquid ~ UL, resulting 
in Pi ~ piU 2 L/Hd- Since the liquid drop will be de- 
formed when P g ~ Pi, one finally obtains the scaling for 
the dimple height and the bubble volume: 



P g ~ LrjgU r / H\, one immediately obtains 



H d ~ RSt 2 ' 3 , 



V h ~ L 2 H d ~ R 3 St 4 / 3 . (5) 



This describes the air bubble in the inertial regime, i.e. 
large impact velocities, in agreement with our experimen- 
tal and numerical findings. 

For large St, corresponding to small impact velocity 
and small droplet radius, capillarity will take over and 
try to smoothen the dimple out. Then the right hand 
side of the Stokes equation ((T|) must be balanced with the 
Laplace pressure 7/c, where k ~ Hd/L 2 is the curvature 
of the dimple. Using once more that the gas pressure 
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as scaling in the capillary regime. Again, this agrees 
well with the experimental and numerical findings. The 
crossover between the regimes, corresponding to the max- 
imal air bubble entrainment, occurs at 



St a ~ Ca 3 J 4 or U 
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Using prefactors obtained from our experimental data in 
figure [H for an ethanol droplet of 0.9 mm radius, this 
translates to an impact velocity U of 0.25 m/s. What is 
the physical reason for the maximum? For higher veloci- 
ties inertia dominates and flattens the droplet at impact. 
For lower velocities and/or smaller droplets the capil- 
lary forces try to keep the drop spherical. In between 
these two regimes the maximal air entrainment under 
the droplet is achieved. 

For many applications air entrainment is undesirable 
and maximal wetting must be achieved. This holds for 
immersion lithography, wafer drying, glueing, agricul- 
tural applications [13, Ell ■ Intriguingly, for inkjet drops 
of radius R ~ 10/im, the optimal velocity according to ([7} 
is approximately 1 m/s. This lies exactly in the range at 
which inkjet usually operates (typically a few m/s), and 
relatively large bubbles will thus be entrapped [lj. For 
immersion lithography the entrapment of even micron- 
sized bubbles can cause practical limitations [l?], fl8| . 
This technology is based on optical imaging of nanoscale 
structures, for which the optics is immersed in water to 
push the limits of spatial resolution. Clearly, it is crucial 
to avoid bubbles or to minimize their size, which also 
has bearing in cleaning and drying of wafers. Ideally, 
one should stay as far as possible from the optimal air 
entrainment impact velocity. Our findings will help to 
achieve this goal and thus optimal wetting. 
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FIG. 2: Maximum entrapment of air bubbles, (a) Dimple height Hd and (b) entrained bubble volume Vb as functions of the 
impact velocity U (upper axes) and Stokes number St (lower axes). The shape of the air layer can be characterized by the 
dimple height Hd and the lateral extension L. Blue circles correspond to high-speed color interferometry measurements, red 
squares correspond to numerical simulations. The straight lines correspond to the derived scaling laws in the capillary regime 
(solid) and inertial regime (dashed) with the respective scaling exponents. The fitted prefactors for the scaling law of the 
dimple height are order unity, 0.3 and 3 respectively for the capillary and inertial regimes. The prefactors for the scaling laws 
of the volume, being the third power of length scales, are larger, namely 9 and 170, respectively. 
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FIG. 3: Comparison of experimental (blue) and numeri- 
cal (red) dimple profiles for two different impact velocities; 
£/=0.2 m/s (St=1.28- 10 -4 ; crossover regime) and C/=0.7 m/s 
(St=3.66 • 10 -5 ; inertial regime). 
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